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ABSTRACT 

There exist sound literature and algorithms for computing 
Liouvillian solutions for the important problem of linear 
ODEs with rational coefficients. Taking as sample the 363 
second order equations of that type found in Kamke's book, 
for instance, 51% of them admit Liouvillian solutions and so 
are solvable using Kovacic's algorithm. On the other hand, 
special function solutions not admitting Liouvillian form ap- 
pear frequently in mathematical physics, but there are not 
so general algorithms for computing them. In this paper 
we present an algorithm for computing special function so- 
lutions which can be expressed using the 2F1, 1F1 or 0F1 
hypergeometric functions. The algorithm is easy to imple- 
ment in the framework of a computer algebra system and 
systematically solves 91% of the 363 Kamke's linear ODE 
examples mentioned. 

Categories and Subject Descriptors 

I.l [Symbolic and algebraic manipulation] : Algorithms. 

General Terms 

Algorithms, design, theory. 

Keywords 

Linear ordinary differential equations, Non-Liouvillian solu- 
tions, hypergeometric solutions. 

Introduction 

Given a second order linear ODE 

y" + A(x)y' + B(x)y = (1) 

where the quantity 1 A' /2 + A 2 /A — B is a rational function of 
x, the problem under consideration is that of systematically 

1 This quantity is an inva rian t under transformations of the 
dependent variable - see 11211 . 
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computing solutions for this ODE even when the solutions 
admit no Liouvillian form 2 . 

The first thing to note is that non-Liouvillian solutions 
which are representable symbolically not as unknown infi- 
nite sums can be represented using special functions, e.g. 
Bessel, Hermite or Legendre functions 1 . In turn, these 
and most of the special functions frequently appearing in 
mathematical physics happen to be particular cases of the 
pF q hypergeometric function for p equal to 0, 1 or 2 and q 
equal to or 1 (see 0). For example, the Bessel functions 
can be expressed in terms of 0F1, all cylindrical functions 
as well as the Hermite, Laguerre, Whittaker and error fam- 
ily of functions can be expressed in terms of 1F1, and all 
Chebyshev, Gegenbauer, Jacobi, Legendre and some others 
can be expressed in terms of 2F1. 

One natural approach is then to directly attempt the com- 
putation of hypergeometric function solutions of these 0F1, 
1F1 and 2F1 types, since in this way we cover at once so- 
lutions involving all the related special functions. Such an 
approach was developed during the year 2001 (see 0), it 
became the main algorithm of the Maple computer algebra 
system for this type of problem since then and it is the sub- 
ject of this paper. The algorithm consists of an equivalence 
approach to the p F q differential equations, is formulated in 
sec.Q 2 and 3, and computes solutions of the form 

v (--;••; fSrf) ( 2 ) 

where P(x) is an arbitrary function and {a, f3, 7, 5, k} are 
constants. 

ft is important to note that the idea of seeking hyperge- 
ometric function solutions for Q or using an equivalence 
approach for that is not new. In '89 Kamran and Olver 
|lj showed how to use an equivalence approach to compute 
Bessel function solutions to eigenvalue problems. Hyper- 
geometric solutions were also discussed by Petkovsek and 
Salvy in '93. Some of the more recent developments 
were presented as computer algebra algorithms too. For 
instance, a classic invariant theory approach was presented 
during 2000 by von Biilow in [5]; in 2001 Willis \7\ presented 
a semi-heuristic algorithm for computing special functions 

2 Functions that can be expressed in terms of exponen- 
tials, integrals, and algebraic functions, are called Liouvil- 
lian functions. The typical example is exp( J R(x), dx) where 
R(x) is rational or an algebraic function representing the 
roots of a polynomial. 



solutions. In 2002 Bronstein and Lafaille [H] presented an 
approach for resolving an equivalence under rational trans- 
formations, between two linear equations in normal form, 
whenever one of them has an irregular singularity 3 . 

There is natural intersection between what these algo- 
rithms can solve but none can claim to extensively cover 
the portions of the problem covered by all the others. If 
compared with the algorithm presented in this paper - we 
called it hyper3 - these other algorithms, both those devel- 
oped before and after hyper3: 

• Do not resolve in a systematic manner all of the 2F1, 
1F1 and 0F1 equivalences; 

• Do not handle the problem of an invariant involving 
fractional or abstract powers; 

• Do not explore automorphisms to avoid uncomputed 
integrals in the solution. 

Also, hyper3 does not require solving systems of algebraic 
equations nor computing Groebner basis nor running differ- 
ential elimination processes nor eliminating parameters by 
composing resultants (all of them expensive computational 
processes), thus resulting in a fast and smooth algorithm 
with little computational cost. These facts, combined with 
the range of problems it solves, for instance taking Kamke's 
book |12| as a testing arena, are at the base of the role hy- 
per3 has today in the Maple differential equation libraries. 

1. COMPUTING 2 Fi, iFx AND Fi HYPER- 
GEOMETRIC SOLUTIONS 

To compute p F q solutions to the idea is to formulate 
an equivalence approach to the p F q underlying hypergeo- 
metric differential equations; that is, to determine whether 
a given linear ODE can be obtained from one of the 2F1, 
1F1 or 0F1 ODEs, respectively given by 



x) y" + ((a + b+ 1) x - c) y + bay = 0, 
xy" + (c-x)y' -ay = 0, 
xy +cy -y = 0, 



(3) 



where {a, b, c} are arbitrary constants, by means of a trans- 
formation of a certain type. If so, the solution to the given 
linear ODE is obtained by applying the same transformation 
to the solution of the corresponding p F q ODE above. 

This approach of course also requires determining the val- 
ues of the hypergeometric parameters {a, b, c} for which the 
equivalence exists, and it is clear that its chances of success 
depend crucially on how general is the class of transforma- 
tions being considered. For instance, one can verify that for 
linear transformations 4 



x-*F(x), y~*P(x)y 



(4) 



with arbitrary F(x), P(x), the problem is too general in that 
to solve it requires solving first the given ODE, so that the 
approach is of no practical use 0. 

3 That also leads to 1F1 solutions of the form including 
its particular 0F1 case, whenever the point of application of 
1F1 is rational in the independent variable. 
4 The problem of equivalence under transformations {x — > 
F(x), y — > P(x) y + Q(x)\ for linear ODEs can always be 
mapped into one with Q(x) — 0, see [Hj- 



The transformations considered in this work are 



ax k +/3 
7 x K + 



(5) 



with P(x) arbitrary and {a, f3, 7, 8, k} constant with respect 
to x. These transformations, which do not conform a class in 
the strict sense 5 , can be obtained by sequentially composing 
three different transformations each of which does constitute 
a class. The sequence starts with linear fractional - also 
called Mobius - transformations 



M 



ax + /3 



7 x + S 

is followed by power transformations 



(6) 



x->x k , (7) 

and ends with linear homogeneous transformations of the 
dependent variable 



V 



Py 



(8) 



So, we are talking of an algorithm that systematically com- 
putes, when they exist, solutions of the form 



ax k +p 
7 x k + S 



(9) 



where p F q is any of 2F1, 1F1 or 0F1. 



p(x) y of the depen- 



1.1 Transformations ; 
dent variable 

The first thing to note is that transformations of the form 
© can easily be factored out of the problem: if two equa- 
tions of the form Q, with coefficients {A(x),B(x)} and 
{C(x), D(x)} respectively, can be obtained from each other 
by means of ©, the transformation relating them is com- 
putable from these coefficients. For that purpose, we rewrite 
both equations in normal form, for instance for Q use 



y — ue 



■ JA/2dx 



to obtain 



A' A 2 

1 — B ]u 

2 4 



(10) 



(11) 



and the transformation relating the two hypothetical ODEs 
exists when the two normalized equations are equal; the 



transformation relating them being y 
what follows we will refer to 



J(C-A)/2dx 



111 



Kx) 



A' A 2 



B, 



(12) 



the coefficient of u in l|ll|l . as the invariant |10j . regardless of 
the fact that this object is only an absolute invariant under 
® and not under P or 0. 



5 By class of transformations we mean a set of transforma- 
tions closed under composition. 



1.2 Transformations x -> f(x) of the indepen- 
dent variable 

By changing x — > F(x) in Q, the invariant Ii of the 
changed ODE can be expressed in terms of the invariant 7o 
of JTJ by 



(13) 



(14) 



h(x) =F' 2 I (F(x)) + S(F') 
where S(x) is the Schwarzian |11) 

3F" 2 F'" 

The form of S(F') is particularly simple when F(x) is a 
power transformation (see (|2SJ) an d also when F(x) is a 
Mobius transformation (JSJ, in which case S(F') — 0. These 
are key facts permitting a simple formulation and resolution 
of the equivalence. 

2. MOBIUS TRANSFORMATIONS AND A 
CLASSIFICATION OF SINGULARITIES 

The first ODE in ^ has 3 regular singularities, at 0, 1 
and oo. The second ODE in JHJ, also known as the confluent 
hypergeometric equation, has a regular singularity at and 
an irregular one at oo. The third ODE in © also has one 
regular and one irregular singularity at and oo, but we 
considered the case separately in order to obtain solutions 
directly expressed in terms of simpler (Bessel) functions. As 
we shall see, the structure of the singularities of these equa- 
tions is a key for resolving related equivalences and Mobius 
transformations preserve that structure. These transforma- 
tions only move the location of the poles. For example, the 
oFi hypergeometric equation 



where all {cji, tjj} can be expressed in terms of {a, b, c} and 
{a, 13, 8, 7} respectively entering the 2F1 equation © and 
the transformation The invariant of the transformed 

1F1 equation has the form 



7 rFr = 



0J2X +2luix + luq 



(osx + a-4) 2 (a 5 x + CTe) 4 
and that of the transformed 0F1 equation has the form 



(18) 



7 oFi = 



LUlX + LOQ 



(03X + <t 4 ) 2 (a 5 x + cto)' 5 



These transformed invariants are all of the form 



UZi (a-ix + bi) 



(19) 



(20) 



Cancellations between factors in the numerator and denom- 
inators of 12(11 may also happen and, independent of that, 
some coefficients {ai,Ci} can be zero 6 . So the degrees with 
respect to x of the numerators and denominators of 117L 
1181 and 11911 can be lower than the maximum implicit by 
these equations; in this way the problem splits into cases. 

Taking these possible cancellations into account, from the 
structure of the invariants 1171 . 118H and 119L the different 
cases for each of the 2F1, 1F1, 0F1 classes were determined. 
With this classification in hands, from the knowledge of the 
degrees with respect to x of the numerator and denominator 
of the invariant I20H of a given ODE, one can tell whether or 
not it can be obtained from the 2F1, 1F1 or 0F1 equations 
© using These observations can be summarized in a 
classification table as follows, using the symbol 



xy" + cy'-y = (15) 

has one regular singularity at the origin and one irregular 
at infinity. The transformed ODE, obtained from 115H by 
means of @ 



y n + (g (Sc + 2yx)+y (2 - c) (3) 



(ax + /3)(~/x + S) 

(a5- 1 f3) 2 
(yx + S) 3 (ax + f3) 



(16) 



also has one regular and one irregular singularity, respec- 
tively located at —(3/ a and —5/^. In the case of the 2F1 
equation (see ©), under the three regular singularities 
move from {0, 1, 00} to {—S/"f, —(3/a, (8 — (3)/ (a — 7)}. So, 
from the structure of the singularities of an ODE, not only 
one can tell with respect to which of the three differential 
equations Q could the equivalence under @ be resolved, 
but also one can extract information regarding the values of 
the parameters {a, (3, 7, 8} entering the transformation. 

Reversing the line of reasoning, through Mobius transfor- 
mations one can formulate a classification of singularities 
of the linear ODEs "equivalent" to the p F q equations JjJJ, 
based on how the invariant of each of these equations is 
transformed. Concretely, after transforming the 2F1 equa- 
tion, the invariant of the resulting equation has the form 



[< p, [qi*,q 2 *,...,q„*]] 

where p is the degree in x of the numerator of 1201 1 and qi are 
the powers of the factors entering the denominator of it. The 
symbol <, when present, refers to the value of p (can be less 
or equal to). The symbol *, when present, means there can 
be factors canceling between numerator and denominator, so 
that the actual value of the related qi can be lower (provided 
p is also lower by the same amount). For example, 



[< 2*, [2*, 2*]] 



(21) 



represents the following possible seven different "lists of val- 
ues" (herein referred as cases) for the degrees of the numer- 
ator and denominator of the invariant 

[2*, [2*, 2*]] = [2,[2,2]],[1,[1,2]],[0,[1,1]],[0,[0,2]] 
[1*,[2*,2*]] = [1,[2,2]],[0,[1,2]] 
[0,[2,2]] 

(22) 

With this notation, the classification of all the possible cases 
equivalent to the 2F1, 1F1 and 0F1 equations under Mobius 
transformations is as shown in Table 1. 



7 2 Fi = 



L02X + 2ui\X + LO(j 



{o\X + 02) ((T3X + 04) (0-5X + 06) 



(17) 



6 Provided that, in (|SJ, a8 — 7/3 7^ and also that in Q the 
invariant remains finite, i.e. its denominator is not zero. 



Class 


Cases 


Number of cases 


2F1 


[<= 2*, [2*, 2*, 2*]], [<= 2*, [2*, 2*]] 


14 


1F1 


[2*, [2*, 4jj, [<= 2, [6jj, [<= 2, [4j|, [2*, [2*]], [2, [0]] 


13 


0F1 


[1*, [2*, 3]], [<= 1, [5]], [<= 1, [3]], [1*, [2*]], [1, [0]] 


9 



Table 1. Classification of linear ODEs equivalent to v F q ODEs under Mobius 



TRANSFORMATIONS x - x* 
INDEPENDENT VARIABLE 



OF THE 



Using the results of the previous sections it is possible to 
resolve the equivalence of a given linear ODE Q and the 
hypergeometric equations Q under compositions of trans- 
formations JHJ of the dependent variable y(x) and Mobius 
transformations @ of the independent variable x. In this 
section a worth additional level of generalization is obtained 
by composing those two transformations with transforma- 
tions x — > x k of the independent variable. 

The first thing to note regarding power transformations 
is that, unlike Mobius transformations, they do not preserve 
the structure of singularities. The change in the invariant 
due to x — > x k , however, has a simple and tractable struc- 
ture. The Schwarzian l|14|l is given by: 



S(F') = 



4x 2 



(23) 



So, the changed invariant I\ shown in 11311 can be expressed 
in terms of 7o by 



x 2 h{x) + \ 
4 



r ! )'l (x k ) + ±)k 2 (24) 



This naturally suggests the introduction of a "shifted" in- 
variant J{x) 



2 1 
Ji{x) = x Ii(x) + - 



(25) 



for which the transformation rule under x — > x has the 
simple form 



Ji(x) — k Jq(x ) 



(26) 



The equivalence of two linear ODEs A and B under x — > 
x k can then be formulated as follows: Given Jia(x) and 
Jib(x), compute Ica and ks entering 126H such that the de- 
grees with respect to x of Jqa{x) and Job{x) are minimized. 
This approach is systematic: equations A and B are related 
through power transformations only when Joa ~ Job and, 
if so, the mapping relating A and B is just x — > x kA ~ kB . 

The computation of k minimizing the degrees of Jo in 1261 
is formulated as follows. Given the set 



i — 1 to m 



(27) 



of (possibly rational) numbers entering as exponents in the 
powers of the independent variable found in J%, compute the 
smallest rational number k such that multiplying by it each 
element of A, all of them become integers. Then the value 
of k minimizing the degrees of Jo is k = 1/k. 



4. SUMMARY OF HYPER3 - EXAMPLES 

An itemized description of the algorithm, discussed in the 
previous subsections to resolve the equivalence proposed in 
the introduction, is as follows. 

1. Rewrite the given equation Q we want to solve in 
normal form 



y" = K x )v 

where I(x) is the invariant 11211 . 



(28) 



2. Compute Ji(a;), the shifted invariant I25L and use 
transformations x — > x k to reduce to the integer mini- 
mal values the exponents of powers entering Jo(a;); i.e., 
compute k and with it compute Jo (x) in 126H . 

3. From 1251 . compute Io(x) and classify its structure of 
singularities according to Table 1, to tell whether an 
equivalence under Mobius transformations is possible 
and to which of the 2F1, 1F1 or 0F1 equations 

4. When the equivalence is possible, from the singular- 
ities of Io{x) and by comparing it with the invariant 
121 H of the transformed p F q equation 7 , compute the 
parameters {a, 6, c} entering the p F q equation such 
that the equivalence exists as well as the parameters 
{a, j3, 7,(5} entering the Mobius transformation @. 

5. Compose the three transformations to obtain one of 
the form 



ax k + P 
~yx k + S ' 



y -» P(x)y 



mapping the p F q equation involved into the ODE be- 
ing solved. 

6. Apply this transformation to the known solution of the 
p F q equation resulting in the desired ODE solution. 

An example of the 2F1 class 

Consider the second order linear ODE 



V = 



2(i/-(j)i 2 -3i 4 -2(ji + v)-1 



v (i/ + 2fc+l)) 



(29) 



This equation has regular singularities at {0, 1, — 1, i, — i}. 
Following the steps outlined in the Summary, we rewrite 
the equation in normal form and then compute the value 

7 At this point, Jo(x) and the shifted invariant of the 
p F q equation have the same degrees. 



of k leading to an equation with minimal degrees for the 
powers entering Jo (x) in P^ . The value found is k = 2. So, 
using 8 



t = i 2 , w = VSeV 2 <— 5 ) J y, (30) 

the given equation Ij29^ can be obtained from 



in a Pentium IV, 2 GigaHertz computer. The Maple com- 
mand line to compute this solution directly using hyper3 is: 
> dsolve (ode , [hyper3] ) ; . 

An example of the 1F1 class 

As an example which also requires an extension of the 
algorithm to handle symbolic powers in the invariant 112L 
consider Kamke's second order linear equation 2.15: 



„ (m 2 + 2 {v 2 - ^ - 2)) t 2 + 2 U 2 - v 2 ) t + n (fi + 2) 

II — — 11 

4t 2 (t-l) 2 (t + l) 2 

(31) 

which is in normal form and has an invariant with "minimal 
degrees" with respect to power transformations 0. 

In step 3, analyzing the invariant of I31H (coefficient of u 
in its right-hand-side), the equation has now three regular 
singular points, at {0, 1, —1}. Using the notation of sec . 11. 21 
the degrees with respect to t of the numerator and of each of 
the linear factors entering the denominator are [2, [2,2,2]]. 
The equation matches the classification Table 1 presented 
in sec. ll.2l and is identified as equivalent to the 2F1 equation 
under Mobius transformations @. 

So we proceed with step 4, equating the invariant of 13H 
with the invariant 1171 written in terms of {a, b, c, a, j3, 7, 8}, 
from where we compute the values of the hypergeometric 
parameters {a, b, c} entering the 2F1 equation such that 
the equivalence under Mobius exists, as well as the Mobius 
transformation itself, obtaining 



{a 



b = 



= -fj,} M :=x = 



It 



The transformation mapping the 2F1 equation © at these 
values of the parameters {o, b, c} into l|31|l is then obtained 
composing the Mobius transformation above with one of the 
form computed as explained in sec. ll.ll resulting in 



2t 
FT' 



t^ 2 (t - 



<(„- M -l)/2 



(t + l) 



(H-l)/2 



■ u(t) (32) 



At this point, we have the transformation I32H mapping 
into 1311 1. and the transformation 13011 . mapping 1311 into 
the equation 1291 we want to solve. Composing these trans- 
formations, in step six we obtain the solution of 1291 



V = 



+ 



(x 2 - 1)- 



V V 

2' 2 



2x 2 



Ch 



(x 2 - l) 1+ f+^ 2 

V 

2 



(33) 



2F1 I ~ + i, i + m+tt; 2 + m; 



2x 2 



1 



where Ci and C2 are arbitrary constants. 

As mentioned in the introduction, an implementation of 
the algorithm being presented is at the core of the current 
Maple ability to solve this type of problem. The time con- 
sumed by this Maple implementation to compute the solu- 
tion 1331 1 performing all the steps mentioned is 0.4 seconds 



8 This transformation is the c omp osition o f t = x k = x 2 with 
a transformation of the form 11 OH so that 1311 is normalized. 



y" + {nx 2 ' T + vx"- 1 )y = Q (34) 

where fi, v and a are constants with respect to x. This 
equation is already in normal form and the shifted invariant 
for it is 



Ji{x) = 1/4 — x 2 (fj, x 2 " + v x a 1 ) 



(35) 

To compute the values of k entering 12611 and leading to 
Jo (a;) with minimized integer powers, in 12711 . instead of 
restricting k to be a rational number, we allow it to de- 
pend on symbolic variables. So we compute k such that 
the set of exponents entering 135L A := {2a + 2, a + 1}, 
becomes a set of integers after multiplying each element of 
it by k, resulting in 9 k = l/(cr + 1). In summary, using 



= x rT+1 , u(t) = x a l 2 2/(2;) j , Kamke's equation 13411 can 

be obtained from the following equation, which is already 
in normal form and has an invariant with minimized integer 
degrees, free of symbolic powers 

(4 M t 2 + 4«/t + g 2 + 2a) 

it = — 5 -u (36) 

4(<r + l) 2 * 2 ^ ; 

Proceeding with step 3, the invariant is the coefficient 
of u in the above and the degrees with respect to t of its 
numerator and factors in its denominator match the Table 1 
of sec. ll.2l identifying 1361 as equivalent to the 1F1 equation 
under Mobius transformations @. 

As in the previous example, in step 4, comparing the in- 
variant of 136H with the invariant 1181 of the transformed 
1F1 equation, we compute the values of the parameters en- 
tering the 1F1 equation © such that the equivalence exists, 
as well as the parameters entering the Mobius transforma- 
tion. Composing all the transformations, we arrive at the 
solution for Kamke's example 2.15 



V = 



a + 1 



1F1 



2^7i(a + l) ; a + 1' a + 1 



x 



(37) 



/ (a + 2) + iv a + 2 2 i ^ a+1 
+ 1 1 ^ 2^7i(a + l) 'a + 1' a + 1 X 



xC 2 



where Ci and C2 are arbitrary constants. The time con- 
sumed by the implementation in Maple to perform these 
steps and return the solution above is again 0.4 seconds, as 
in the previous example. This also illustrates that, for typ- 
ical problems, the additional handling of symbolic powers 
does not imply on any important performance cost. 

9 To perform this computation, it suffices to sequentially take 
the gcd between each of the elements of A. 



5. ON THE COMPUTATION OF THE SEC- 
OND INDEPENDENT SOLUTION 

The algorithm presented is based on computing a trans- 
formation mapping a p F q equation into a given linear ODE, 
then applying that transformation to the solution of the 
pF q equation to obtain the solution for the given problem. 
This process has a subtlety: depending on the values of the 
hypergeometric parameters, we may have only one indepen- 
dent solution available for the p F q equation. In these cases, 
the second independent solution can be obtained through 
integration: if y = S(x) is a solution of then 

r MA(x)dx) 

V = J 2 dxS(x) (38) 

is a second independent solution directly computable from 
S(x) and A{x). 

This approach, however, frequently introduces uncompu- 
table integrals, thus complicating further manipulations and 
undermining the usefulness of the result. As an example of 
this situation, for the 2F1 equation, 



from Table 2 by respectively changing {0, 1, 00} by {A, /x, k}. 
Hence, the solution 140H can be written in different manners, 
by changing the application point of the 2F1 function using 
the Qi, permuting accordingly the parameters {\,[i,k} en- 
tering the 2F1 function and multiplying the result by the 
proper non-constant factor 10 . 

For example, when c is an integer but a + b is not an 
integer, applying g 2 and permuting the parameters fj, <-> A, 
the power x x ~ c entering 1401 1 becomes a power with non- 
integer exponent. Using this mechanism, for at c = 1, 
instead of the solution with integrals 1411 we obtain two 
independent solutions free of uncomputed integrals: 

y = 2 F 1 (a,b;a + b;l-x)C 1 (43) 
+ {x- i) 1 - b ~ a 2 Fr (1 - b, 1 - a; 2 - 6 - a; 1 - x) C 2 

When c and a + b are both integers, #2 does not resolve 
the problem, but if a — b is not an integer then g$ does, 
since it permutes the integer A = 1 — c with the non-integer 
k — a — b. For example, for a — 2/3, b — 1/3, c = 1, 1391 
becomes 



(x 2 - x) y" + ((a + b + 1) x - c) y' + bay = 0, (39) 
the two independent solutions are: 

y = 2F1 (a,b; c; x) Ci (40) 
+ x 1 ~ c 2 Fi (6 - c+ 1, a - c+ 1; 2 - c; x) C 2 

but for c = 1 these two solutions are equal. Using the inte- 
gration recipe 138H . a second independent solution is 



V 



(a + b+ l)x - 1 



dx 



■dx 2F1 (a,b; 1; a;) (41) 



2F1 (a, b; 1; x) 2 

Although the inner integral, with rational integrand, is easy 
to compute, the outer integral, with 2F1 (a, b; 1; x) 2 in its 
denominator, is uncomputable in current computer algebra 
systems. 

The approach used in hyper3 to minimize the occurrence 
of uncomputable integrals consists of exploring the group 
of automorphisms of the 2F1 equation in order to make c 
not an integer when that is possible. Recalling, the group 
elements and their action are 



Group element 



9i ■ x 



92 -x 



1-x 



93 ■ x 



~TJx~ 



94 



x -> 1/(1 



gs : x — > (x - I)/ 2 



Action on the plane 



(0 -> 0, 1 -> l,oo 



"TO" 



1,1 



0, 00 



<x, 1 — > 1, 00 



To" 



w 



1,1 



IE 



00, 1 — » 0, 00 



1 



<?6 : x — > x/(x — 1) (0 — > 0, 1 — > 00, co — > 1) 

Table 2. Group of automorphisms of the 2^1 equation 

These transformations, known to act as permutations on 
the set {0, 1, 00}, also act as permutations on a set {A, fi, ft} 
related to the hypergeometric parameters {a, b, c} by 



A = 1 — c, fj, — a + b 



(42) 



These three parameters are the exponent differences of the 
normal form of the 2F1 equation Q, at {0,1, 00} respec- 
tively. The action of each t/; on these parameters is obtained 



2 y/9 + (2 x - 1) y + (x 2 - x) y" 







(44) 



Applying (73 and permuting the parameters A and k, we ob- 
tain the following two independent solutions free of integrals 



V = 



- 1/3 2 Fi(l/3,l/3; 2/3; 1/x) d 



(45) 



+ x 2/ ' i 2 F 1 (2/3,2/3; 4/3; 1/x) C 2 



When all of c, a + b and a — b are integers, these permuta- 
tions are in principle of no use, but still for some cases the 
solution can be represented free of integrals. This is the case 
of Legendre's equation. Recalling the relationship between 
the associated Legendre function of the first kind and the 
hypergeometric 2F1 function , 



LegendreP (a, b, z) 



(46) 



(2 + 1) 



l/2i> 



■Fi (a+ 1, 



1 - b; (1 - z)/2) 



(z-l) 1/2b F{l-b) 



whenever the group elements of Table 2 can map the 2F1 
function solution into one of the form above, one indepen- 
dent solution can be expressed using LegendreP and the sec- 
ond one is obtained from the first one replacing LegendreP 
by the associated function of the second kind LegendreQ. 
For example, for 



y/4 + (2x- l)y'+ (x 2 - x) 



x) y 







(47) 



we have /j, — k = \ = 0, so c = 1 and both a + b and a — b 
are integers. A solution free of integrals is 



y = LegendreP (-1/2, 2 a; - 1) Ci 

+ LegendreQ (-1/2, 2 a; - 1) C 2 



(48) 



10 These multiplicative factors are different for each gi; we 
omit them here for brevity. 

11 We use here the Maple convention for the branch cuts of 
LegendreP; the idea being discussed is independent of that. 



Conclusions 

In this presentation we discussed an algorithm for second 
order linear ODEs, we called it hyper3, for computing non- 
Liouvillian solutions by resolving an equivalence to the 2F1, 
1F1 and 0F1 equations. Taking Kamke's book as testing 
arena, this algorithm is the most successful one of the cur- 
rent set of linear ODE algorithms of the Maple system. From 
the 363 corresponding examples of Kamke's book having ra- 
tional coefficients, hyper3 alone solves 331 (91 %), followed 
by Kovacic's algorithm solving 181 (50 %). Moreover, from 
these 181 examples admitting Liouvillian solutions, hyper3 
solves 163 (90 %). 

The fact that, for 90% of these equations admitting Li- 
ouvillian solutions, the solution can also be computed as a 
hypergeometric one of the form @ is a good indication that 
the restriction used to make the algorithm feasible is appro- 
priate. The fact that around one half of Kamke's examples 
only admit special function solutions of non-Liouvillian form 
also illustrates the relevance of this type of solution in the 
general framework of linear ODE problems popping up in 
applications. 

Despite the simplicity of the approach, till the end of 2001, 
when the routines for this algorithm were developed, no 
equivalent or similar algorithms were available in any of the 
Axiom, Maple, Mathematica, MuPAD or Reduce computer 
algebra systems (CAS). These CAS failed in computing spe- 
cial function solutions but for occasional success, e.g., by 
previous to hyper3 Maple routines able to resolve an equiv- 
alence under only power transformations of the form J7J 
|13|. or an equivalence under only Mobius transformations 
and only with respect to the 2F1 class |14|. 

Since at the core of hyper3 there is the concept of singular- 
ities, two natural extensions of this work consist of applying 
the same ideas to compute solutions for linear ODEs of or- 
der three and higher I15| and for second order equations 
of Heun type. The latter have four regular singular points 
or any combination of singularities derived from that case 
through confluence processes 116) : one example of these are 
Mathieu equations. Related work is in progress |17II18| . 
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